UNCLASSIFIED 


AD  NUMBER 


AD883322 


NEW  LIMITATION  CHANGE 
TO 

Approved  for  public  release,  distribution 
unlimited 


FROM 

Distribution  authorized  to  U.S.  Gov't, 
agencies  only;  Test  and  Evaluation;  APR 
1971.  Other  requests  shall  be  referred  to 
Naval  Ship  Research  and  Development 
Center,  Washington,  DC  20034. 


AUTHORITY 


NSRDC  ltr  dtd  19  Mar  1974 


THIS  PAGE  IS  UNCLASSIFIED 


NAVAL  SHIP  RESEARCH  AND  DEVELOPMENT  CENTER 


co 

CO 

00 

00 


Washington,  D.C.  20034 


Ww 


HYDRODYNAMIC  PERFORMANCE  OF  A  PROPOSED  CRUISING 
DUCTED  PROPELLER  FOR  SUBMARINES 


,  % 


John  L.  Beveridge 


Distribution  limited  to  U.S.  Government  agencies  only; 
Test  and  Evaluation;  28  Nov  1970.  Other  requests  for 
Oris  document  must  be  referred  to  NSRDC,  Code  500. 


April  1971 


SHIP  PERFORMANCE  DEPARTMENT 
RESEARCH  AND  DEVELOPMENT  REPORT 


D  D  C 


Miv  io  197 J 


Report  3374 


r 


The  Naval  Ship  Research  and  Development  Center  is  a  U.S.  Navy  center  for  laboratory  effort 
directed  at  achieving  improved  sea  and- air  vehicles.  It  was  formed  in  March  1967  by  merging  the 
David  Taylor  Model  Basin  at  Carderock,  Maryland  and  the  Marine  Engineering  Laboratory  (now 
Naval  Ship  R  &  D  Laboratory)  at  Annapolis,  Maryland.  The  Mine  Defense  Laboratory  (now  Nava! 
Ship  R  &  D  Laboratory)  Panama  City,  Florida  became  part  of  the  Center  in  November  1967. 

Naval  Ship  Research  and  Development  Center 
Washington,  D.C.  20034 


MAJOR  NSRDC  ORGANIZATIONAL  COMPONENTS 


DEPARTMENT  OF  THE  NAVY 
NAVAL  SHIP  RESEARCH  AND  DEVELOPMENT  CENTER 

WASHINGTON,  D.C.  200  34 


HYDRODYNAMIC  PERFORMANCE  OF  A  PROPOSED  CRUISING, 
DUCTED  PROPELLER  FOR  SUBMARINES 


by 


John  L.  Beveridge 


Distribution  limited  to  U.S.  Government  agencies  only; 
Test  and  Evaluation;  28  Nov  1970.  Other  requests  for 
this  document  must  be  referred  to  NSRDC,  Code  500. 


April  1971 


Report  3374 


TABLE  OF  CONTENTS 


Page 


ABSTRACT  . - .  1 

ADMINISTRATIVE  INFORMATION  . .  1 

INTRODUCTION  .  1 

GEOMETRY  .  2 

HULL  .  2 

DUCT  AND  PROPELLER  . 2 

APPROACH  .  4 

DESIGN  INFORMATION  . 5 

COMPUTER  PROGRAMS  .  5 

PROCEDURE  AND  NUMERICAL  DATA  .  5 

THEORETICAL  PERFORMANCE  . 6 

DUCT  PRESSURE  DISTRIBUTION  . . 6 

EFFICIENCY  OF  DUCT  AND  PROPELLER  SYSTEM  .  7 

PROPULSION  TEST  RESULTS  AND  DISCUSSION  .  8 

CONCLUSIONS  .  11 

ACKNOWLEDGMENTS  .  11 

REFERENCES  .  22 

LIST  OF  FIGURES 

Figure  1  -  NSRDC  Model  4620  with  Ducted  Propeller  4271  .  12 

Figure  2  -  Comparison  of  Hull  Shapes  for  NSRDC  Model  4620  and  for  the 

Point  Source  and  Line  Sink  of  Reference  1  .  13 

Figure  3  -  Cruising  Duct  2.1  13 

Figure  4  -  Duct  Section  and  Dimensions  .  14 

Figure  5  -  Details  of  Propeller  4271  .  14 

Figure  6  -  Hull  Wake  and  Total  Wake  of  Hull  Plus  Duct  .  15 

Figure  7  —  Additional  Angles  of  Attack  for  Blade  Sections  of 

Propeller  4271  16 


ii 


Page 


Figure  8  -  Final  Blade-Pitch  Distribution  for  Propeller  4271  .  17 

Figure  9  -  Pressure  Distribution  on  Duct  .  18 

Figure  10  -  Ducted  Propeller  Efficiencies  versus  RPM  .  19 

Figure  11  -  Model  4620  Overload  and  Underload  Propulsion  with 

Ducted  Propeller  4271  . ... .  20 

Figu'.e  12  -  Power  and  RPM  Curves  for  a  200-Foot  Prototype  of 

Model  4620  .  21 


LIST  OF  TABLES 

Table  1  -  Offsets  and  Particulars  for  NSRDC  Model  4620  .  3 

Table  2  -  Duct  Camber  Ordinates  .  4 

Table  3  -  Comparative  Propulsion  Data  for  a  200-Foot  Prototype  of 

Model  4620  at  V  =  29.3  knots  . . .  9 


NOTATION 


cT 

(CTh)  s 


TP 


Propeller  expanded-blade  area 
Propeller  disk  area 
Pressure  coefficient  p/q 

1  i 

Propeller  power  coefficient  PD  /  —  pAQV  J 

2 

Total  resistance  coefficient  RT  /  — pSP2 

1  2 

1  9 

Propeller  thrust-loading  coefficient  T p  /  —  p  AQV  z 

,  1  , 

Propeller  thrust-power  coefficient  T p  Va  /—  p  AQ  V  3 


c 

D 

J< 

L 

i 


n 

P 

Pr 


V 

Q 


9 

ft 


rpm 


S 

shp 


Duct  chord  length 

Propeller  diameter  or  maximum  diameter  of  a  body  of  revolution 
Propeller  advance  coefficient  V  /  nD 

Body  length 

Propeller  blade  section  length 
Propeller  rate  of  revolution  in  rps 
Propeller  geometric  pitch 
Power  delivered  to  propeller  2ttQti 

Effective  (tow  rope)  power  RT  V 

Local  static  pressure,  excluding  hydrostatic  pressure 
Propeller  torque 

„  1  , 

Free-stream  stagnation  pressure  —  pVz 

2 

Propeller  radius  or  radius  in  general 
Total  resistance 

Propeller  revolutions  per  minute 
Camber  offset,  measured  from  duct  axis 

Wetted  surface 
Shaft  horsepower 


IV 


a.  Wx ffykAi' ^ -'  "■ci>-.;:""“  » .-v ??»  V***JM«y ViW-Wv: »A  ^  ■* */Z4 


T 


Thrust 


£ 

V 

V 


l-wc 

l-wx 

X 

Xc 

X 

Y 

y 

z 


Pi 

( 

^B 

VD 

VH 

T]t 

9 

d 

h 

i 

P 


Thrust-deduction  coefficient 
Ship  speed 

Speed  of  advance  (1  -  w0)P,  includes  hull  and  duct  velocities 
where  appropriate 

Axial  component  of  local  velocity,  includes  hull  and  duct 
velocities  where  appropriate 

Propeller  effective  inflow  velocity  ratio  V a  /  V 
Propeller  local  inflow  velocity  ratio  V %  /  V 

Axial  coordinate,  origin  at  bow  for  hull,  at  leading  edge  for  duct 
Nondimensional  duct  length 

Propeller  radius  fraction  or  nondimensional  length 
Offset  of  meridian  profile  for  body  of  revolution 
Nondimensional  offset  Y  /  D 
Number  of  propeller  blades 

Ideal  angle  of  attack  due  to  loading  in  degrees 

Angle  of  attack  due  to  thickness  in  degrees 
Hydrodynamic  pitch  angle 
Section  ratio  of  drag  to  lift 
Propeller  efficiency  CTp  /  Cp$  ,  Equation  (6) 

Quasi  propulsive  coefficient  Pp  /  PD 

Hull  efficiency  (1  -  t)  /  (1  -  wQ) 

Total  efficiency  of  duct  plus  propeller  VBll  +  (Td)/(Tp)) 

Mass  density 


Subscripts 

Duct  s 

Propeller  hub  Q 

Nonviscous  B 

Propeller 


Ship  speed 

Source 

Sink 


v 


ABSTRACT 


A  theoretical  and  experimental  investigation  was  performed  to  determine 
the  steady-state  propulsion  characteristics  of  a  cruising,  ducted-propeller  system 
for  a  submersible  body  of  revolution.  The  most  important  aspects  of  the  design 
method  and  numerical  procedures  are  discussed.  The  theoretically  predicted  and 
experimental  performance  of  the  ducted  propeller  showed  good  agreement  on  power. 
However  the  experimental  rpm  was  6  percent  higher  than  predicted;  this  difference 
is  attributed  to  the  fact  that  the  propeller  was  underpitched  due  to  an  error  in  input 
to  the  computer  program.  The  subject  propulsion  system  was  found  to  be  less  effi¬ 
cient  than  an  optimum  unshrouded,  wake-adapted  propeller,  but  the  optimum  rpm- 
diameter -relationship  was  not  determined  for  the  ducted  propeller.  (  ^ 


ADMINISTRATIVE  INFORMATION 


The  work  was  authorized  under  the  Naval  Ship  Systems  Command  General  Hydro¬ 
mechanics  Research  Program  and  was  funded  under  Subproject  SR009-0101,  Task  0101. 


INTRODUCTION 

Ducted  propellers  offer  an  alternative  propulsion  system  to  conventional  stern  propel¬ 
lers  for  ships  and  submarines.  Specific  types  of  ducted  propellers  have  been  used,  depending 
on  the  design  conditions  to  be  met.  Kort  nozzle  and  pumpjet  designs  are  well-known  examples 
of  flow-accelerating  and  flow-decelerating  types,  respectively.  The  primary  objective  of  the 
present  investigation  was  to  study  the  feasibility  of  using  a  so-called  cruising,  ducted- 
propeller  propulsion  system  on  a  high-speed  submarine.  Kriebeland  Mendenhall1  have  made  a 
theoretical  analysis  of  the  flow  over  an  underwater  hull  (the  same  hull  form  is  treated  in  this 
report)  with  a  stern-mounted  ducted  propeller  (flow  accelerating  duct).  The  optimum  duct  cam¬ 
ber  (largest  duct  thrust  without  duct  flow  separation)  derived  in  Reference  1  was  used  in  con¬ 
junction  with  the  Morgan2  ducted-propeller  theory  and  a  computer  program  reported  by  Caster3 
to  design  and  predict  the  performance  of  a  cruising,  ducted  propeller. 

A  cruising,  ducted  propeller  can  be  loosely  defined  in  terms  of  the  shape  of  the  duct- 
mean  camber,  that  is,  a  camber  with  the  following  properties: 

1.  No  leading  edge  suction. 

2.  No  hull-duct  interference  (no  “thrust  deduction”)  except  that  due  to  the  propeller. 

S.  Maximum  duct  thrust  at  the  onset  of  separation. 
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References  are  listed  on  page  22. 
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If  a  ducted  propeller  possessing  these  properties  can  be  designed  to  compare  favorably  with 
the  propulsive  efficiency  of  an  unshrouded  stern-mounted  propeller,  then  other  possible  ad¬ 
vantages  can  be  utilized  without  penalty  in  absorbed  power.  It  is  visualized  that  the  stern 
control  surfaces  (rudders  and  diving  planes)  would  be  mounted  externally  on  the  duct.  Possi¬ 
ble  collateral  benefits  include  (1)  a  reduction  in  noise  generation  and  induced  vibration  of 
the  hull  and  (2)  an  added  stabilizing-fin  effect  from  the  duct.  This  report  presents  and  dis¬ 
cusses  the  theoretical  design  techniques  and  the  experimental  results  used  in  a  hydrodynamic 
evaluation  of  an  optimum  cruising,  ducted  propeller  on  NSRDC  Hull  Model  4620.4 

GEOMETRY 

HULL 

NSRDC  Model  4620  is  a  15.0-ft  submersible  hull  (free-flooding)  constructed  of  fiber¬ 
glass  and  equipped  with  standard  NSRDC  instrumentation  for  measuring  body,  duct,  and  pro¬ 
peller  forces.  Table  1  gives  the  offsets  and  principal  dimensions  of  the  model  hull.  Figure 
1  shows  several  views  of  the  model  with  the  ducted  propeller  mounted  in  place  for  testing. 

Reference  l  approximated  the  profile  of  the  body-of-revolution  hull  by  a  point-source, 
line-sink  distribution.  The  strength  of  the  line  sink  was  assumed  to  be  affine  to  the  slope 
of  the  hull  sectional  area  curve.  Since  the  singularities  used  to  generate  the  hull  shape  are 
an  important  part  of  the  analytical  derivation  of  the  duct  camber,  it  is  interesting  to  compare 
the  hull  shape  for  a  point  source  and  line  sink  with  that  of  NSRDC  Model  4620.  Note  in  Fig¬ 
ure  2  that  the  chosen  singularity  system  accurately  approximated  the  actual  hull  shape  in  the 
vicinity  of  the  duct.  The  analysis  of  Reference  1  was  insensitive  to  the  shape  of  the  bow 
portion  of  the  hull. 


DUCT  AND  FROPELLER 


Figure  3  shows  the  aluminum  duct.  An  exploded  view  is  included  to  show  the  various 
components  of  the  ducted  propeller;  major  features  include  propeller,  cruciform  supporting 
struts,  and  the  “backbone  eiement,”  which  contains  propeller-shaft  bearings  and  differential- 
reluctance  block  gages  for  measuring  the  duct  axial  force.  Transparent  plastic  windows  were 
provided  in  the  duct  surface  for  viewing  the  propeller,  should  cavitation  tests  be  performed. 

Duct  camber  ordinates  are  given  in  Table  2.  The  NACA  16-009  thickness  form  was 
added  to  the  camber  to  obtain  the  duct  section.  The  resulting  section  is  shown  in  Figure  4 
together  with  pertinent  duct  dimensions. 

Figure  5  shows  seven-bladed  Propeller  4271  which  was  designed  for  the  duct.  Ex¬ 
panded  blade-section  lengths  were  calculated  from  the  relation 


l/D  = 


Z  (1  -  xh) 
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TABLE  1 

Offsets  and  Particulars  for  NSRDC  Model  4620 


X 

X 

inches 

y 

}' 

inches 

0.00 

000.0 

0.0000 

0.000 

0.02 

3.6 

0.1427 

3.500 

0.04 

7.2 

0.2029 

4.977 

0.06 

10.8 

0.2490 

6.108 

0.08 

14.4 

0.2873 

7.047 

0.10 

18.0 

0.3200 

7.850 

0.12 

21.6 

0.3485 

8.549 

0.14 

25.2 

0.3734 

9.160 

0.16 

28.8 

0.3953 

9.697 

0.18 

32.4 

0.4145 

10.17 

0.20 

36.0 

0.4312 

10.58 

6.22 

39.6 

0.4457 

10.93 

0.24 

43.2 

0.4581 

11.24 

0.26 

46.8 

0.4687 

11.50 

0.28 

50.4 

0.4775 

11.71 

0.30 

54.0 

0.4848 

11.89 

0.32 

57.6 

0.4905 

12.03 

0.34 

61.2 

0.4947 

12.13 

0.36 

64.8 

0.4977 

12.21 

0.38 

68.4 

0.4994 

12.25 

0.40 

72.0 

0.5000 

12.27 

0.42 

75.6 

0.4995 

12.25 

0.44 

79.2 

0.4979 

12.21 

0.46 

82.8 

0.4953 

12.15 

0.48 

86.4 

0.4917 

12.06 

0.50 

90.0 

0.4878 

11.97 

0.52 

93.6 

0.4818 

11.82 

0.54 

97.2 

0.4755 

11.66 

0.56 

100.8 

0.4684 

11.49 

0.58 

104.4 

0.4603 

11.29 

0.60 

108.0 

0.4513 

11.07 

0.62 

111.6 

0.4414 

10.83 

0.64 

115.2 

0.4305 

10.56 

0.66 

118.8 

0.4187 

10.27 

0.68 

122.4 

0.4058 

9.954 

0.70 

126.0 

0.3919 

9.613 

0.72 

129.6 

0.3768 

9.243 

0.74 

133.2 

0.3605 

8.843 

0.76 

136.8 

0.3429 

8.411 

0.78 

140.4 

0.3239 

7.945 

0.80 

144.0 

0.3036 

7.447 

0.82 

147.6 

0.2817 

6.910 

0.84 

151.2 

0.2582 

6.334 

0.86 

154.8 

0.2330 

5.715 

0.88 

158.4 

0.2060 

5.053 

0.90 

162.0 

0.1771 

4.344 

0.92 

165.6 

0.1461 

3.584 

0.94 

169.2 

0.1131 

2.774 

0.96 

172.8 

0.0778 

1.908 

0.98 

176.4 

0.0401 

0.984 

1.00 

180.0 

0.0000 

0.000 

Formula: 

y2  =  a^x  +  a  2  a2  +  +  a^x1  +  +  fig  a:6 

where: 

rjj  =  1.000000 

a2  =  1.137153 

a3  =  -10.774885 
a4  =  +19.784286 
a5  =  -16.792534 
a6  =  +  5.645977 

Wetted  Surface  Coefficient  =  0.7324 
LCB.x  -  0.4456 
L/D  =7.339 

Model  Particulars: 

Length,  ft  =  15.0000 
Diameter,  ft  =  2.044  (24.53  in.) 

Nose  Radios,  ft  =  0.1392  (1.670  in.) 

Tail  Radius,  ft  =  0.0000 
Wetted  Surface,  ft2  =  70.55 
Volume,  ft3  =  29.53 
LCB,  ft  =  6.6840 

Ship  Particulars: 

L  =  200  ft 
77  =  27  ft  3  in. 
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where 


£  is  section  length, 

D  is  propeller  diameter, 

Ag  is  the  expanded  area  of  all  blades, 
A0  is  the  propeller  disk  area, 
x  is  the  radius  fraction, 
xh  is  the  radius  fraction  of  the  hub,  and 
Z  is  the  number  of  propeller  blades. 

The  previous  expression  results  in  a  wide  tip 
propeller  of  the  Kaplan  type.  A  NSRDC 
modified-66-thickness  form,  cambered  with  an 
a  0.8  mean  line,  was  used  for  the  blade 
sections. 


TABLE  2 

Duct  Camber  Ordinates 


x/c 

rs/c 

'0.00 

0.682 

0.05 

0.683 

0.10 

0.679 

0.20 

0.665 

0.30 

0.642 

0.40 

0.614 

0.50 

0.585 

0.60 

0.557 

C.70 

0.532 

0.80 

0.513 

0.90 

0.502 

0.95 

0.500 

1.00 

0.500 

APPROACH 

An  optimum  duct  camber  (Duct  2.1  of  Reference  1)  was  specified.  With  the  duct  cam¬ 
ber  and  duct  size  chosen,  the  problem  was  to  design  and  predict  the  performance  of  the  ducted 
propeller  that  absorbed  minimum  power  at  design  thrust  and  speed.  To  solve  the  problem,  the 
work  of  Morgan  and  Caster2,3,5  was  used.  The  duct  forces,  circulation,  pressure  distribution, 
and  the  effect  of  the  duct  on  the  propeller  were  computed  by  the  linearized  theory  of  the  duct 
presented  by  Morgan.2  Propeller  design  and  performance  were  based  on  the  Lerbs  lifting-line 
theory6  for  moderately  loaded,  finite-bladed  propellers.  In  determining  the  effect  of  the  pro¬ 
peller  on  the  duct,  the  average  axial  and  radial  velocity  components  induced  on  the  duct  by 
the  propeller  were  obtained  from  a  propeller  actuator  disk  theory  as  developed  by  Hough  and 
Ordway.7 

As  mentioned  previously,  the  submersible  hull  form  used  for  the  study  was  identical  to 
the  nodel  hull  for  which  computations  were  performed  in  Reference  1.  Experimental  wake 
data,  needed  in  the  ducted-propeller  design  procedure,  have  been  previously  reported  for  the 
hull  form  used.4  Reference  3  contains  a  complete  discussion  of  the  form  and  nature  of  the 
input  required  for  the  ducted-propeller  computer  program,  the  various  options  available,  the 
output  provided,  and  the  assumptions  of  and  limitations  to  the  theory.  It  seems  desirable  to 
restate  here  the  most  important  assumptions  and  limitations: 

“1.  The  fluid  is  inviscid  and  incompressible  and  no  separation  occurs.  .  .  .  The  viscous 
drag  of  the  duct,  which  includes  both  the  skin-friction  and  pressure  drag,  can  also  be 
calculated  .... 

“2.  The  free-stream  flow  must  be  axisymmetfic. 

“3.  The  duct  is  axisymmetric  and  of  finite  length. 


“4.  The  duct  can  be  represented  mathematically  by  a  distribution  of  ring  vortices  and 
ring  sources  along  a  cylinder  of  constant  diameter.  This  implies  that  the  boundary 
condition  are  linearized.” 

A  parametric  study  was  made  to  optimize  propeller  rpm  for  the  given  duct  diameter. 

The  optimum  rpm-diameter  relationship  would  normally  be  determined  from  a  parametric  study, 
but  the  more  limited  objective  of  evaluating  the  subject  configuration  was  pursued  in  the 
present  case.  Although  the  Lerbs  lifting-line  propeller  theory  was  used  to  design  (in  the 
sense  of  obtaining  the  best  blade  radial-loading  distribution)  and  to  predict  the  performance 
of  a  series  of  propellers,  the  final  propeller  design  for  optimum  rpm  was  based  on  propeller 
lifting-surface  theory  (blade  thickness  included). 

DESIGN  INFORMATION 


COMPUTER  PROGRAMS 

Three  IBM-7090  computer  programs  were  used  to  assist  in  the  design,  performance 
prediction,  and  analysis  of  the  ducted-propeller  propulsion  system. 

Program  1  was  a  composite  program  3  with  several  options  which  could  be  used  either 
to  predict  the  hydrodynamic  characteristics  of  an  axisymmetric  duct  (e.g.,  duct  thrust  and 
duct-pressure  distribution)  or  to  iterate  a  complete  ducted-propeller  design.  Particularly 
useful  were  the  options  which  permitted  including  input  for  an  arbitrary  axisymmetric  veloc¬ 
ity,  including  propeller  (or  otherwise)  induced  velocity  at  the  duct  surface,  obtaining  the 
ideal  angle  of  attack  of  the  duct  section,  and  obtaining  the  velocity  field  inside  the  duct. 

Program  2  was  a  lifting-line  propeller  program8  which  was  used  (1)  to  find  either  the 
optimum  or  nonoptimum  radial  distribution  of  hydrodynamic  pitch  angle  (to  produce  the  de¬ 
sign  thrust  or  power)  and  (2)  to  predict  propeller  performance.  The  program  provided  the 
necessary  input  for  the  final  propeller  design  by  lifting-surface  theory. 

Program  3  was  a  composite  design  method9, 10  which  was  used  to  obtain  the  final 
pitch  distribution  and  camber  ratios  for  the  propeller,  based  on  lifting-surface  theory  with 
blade  thickness  included. 

The  following  sections  discuss  how  these  various  programs  and  options  were  used  in 
the  present  analysis. 

PROCEDURE  AND  NUMERICAL  DATA 

A  thrust-loading  coefficient  C Ths  =  T /%  pAQ  V2  =  0.40  was  used  as  specified  in 

Reference  1  for  design  of  the  ducted  system  (duct  loading  and  propeller  loading  relate  to 
model  scale).  Duct  boundary-layer  characteristics  and  shear  drag  as  well  as  hull  resistance 
and  wake  data  were  calculated  for  the  model  range  of  Reynolds  number.1,4 

With  the  duct  geometry  as  input,  Program  1  was  used  to  obtain  duct  thrust,  duct  pres¬ 
sure  distribution  (discussed  in  the  next  section),  and  the  axial  component  of  the  velocity 
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induced  by  the  duct  at  the  propeller  plane.  Figure  6  shows  curves  of  the  experimental  wake4 
for  Model  4620  hull  and  the  total  propeller  inflow  (hull  plus  duct).  The  option  of  Program  1 
which  iterated  a  complete  ducted-propeller  design  was  used,  in  conjunction-with  Program  2 
for  a  parametric  study  to  optimize  propeller  rpm  for  the  given  diameter  and  thrust  loading. 

The  final  propeller  design  (Figure  5)  was  then  completed  by  using  Program  3. 

Some  important  aspects  of  the  propeller-design  method  according  to  lifting  surface 
theory  will  be  briefly  described.  Essentially  the  study  utilized  a  combination  of  work  by 
Kerw-in  and  Leopold,10  Pien,11  and  Cheng.9  The  chordwise  distribution  of  camber  ratio  was 
obtained  for  nine  radii  by  the  Pien  mathematical  model,  which  used  a  continuous  distribution 
of  vortices.  The  final  pitch  angle  of  each  blade  section  was  obtaine'd  by  adding  two  angles 
of  attack  to  the  hydrodynamic  pitch  angle,  the  ideal  angle  of  attack  due  to  loading  from  lifting- 
surface  theory,  and  an  angle  of  attack  due  to  blade  thickness.  The  additional  angles  are 
show'n  in  Figure  7,  and  the  final  pitch  distribution  is  given  in  Figure  8. 

THEORETICAL  PERFORMANCE 
DUCT  PRESSURE  DISTRIBUTION 

Figure  9  is  a  composite  graph  which  shows  the  variation  of  the  theoretical  pressure 
coefficient  versus  chordwise  position  xc  for  points  on  the  duct  outside  and  inside  surfaces. 
The  curve  shows  the  pressure  distribution  calculated  by  the  present  method  for  the  final 
ducted-propeller  configuration  at  optimum  propeller  rpm.  Data  points  are  shown  for  other  con¬ 
ditions  as  follows: 

1.  Computed  results  given  in  Reference  1,  based  on  a  duct  of  zero  thickness  (indicated 
byp  and  o'). 

2.  Results  computed  by  the  present  method  for  the  final  configuration  but  assuming  a  duct 
of  zero  thickness* (indicated  by  O  and  Cf). 

Results  of  the  two  approaches  to  a  solution  for  the  duct  pressure  distribution  were  in 
fairly  good  agreement  for  the  last  half  of  the  duct  length,  particularly  on  the  duct  outside  sur¬ 
face.  As  expected,  the  pressures  calculated  by  the  present  method  were  usually  closer,  to 
the  data  of  Reference  1  when  duct  thickness  was  zero.  Integration  of  the  present  pressure 
distribution  for  zero  duct  thickness  yielded  a  value  of  0.06  for  T d/T p.  This  compares  to  the 
value  of  0.07  obtained  in  Reference  1,  an  insignificant  difference.  If  thickness  is  included, 
Td/Tp  =  0.04. 

An  increasing  pressure  gradient  with  respect  to  xc  was  indicated  in  the  forebody  region 
of  the  duct  outside  surface.  This  could  lead  to  laminar  separation  with  a  consequent  deteri¬ 
oration  in  duct  performance;  the  increase  in  actual  pressure  drag  would  undoubtedly  offset  any 
reduction  in  the  duct  shear  force.  According  to  Reference  1,  the  transition  from  laminar  to 
turbulent  flow  occurred  at  xc  “  0.15  at  design  condition  on  both  the  inside  and  outside  duct 
surface.  Theoretically,  at  design  condition  the  flow  is  not  separated  from  either  the  inside 
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or  outside  duct  surface  for  the  subject  duct.  As  will  be  discussed  later,  the  propulsion  test 
included  a  large  range  of  loadings,  and  the  resulting  forces  undoubtedly  reflected  a  flow  that 
separated  far  forward  of  the  duct  trailing  edge  (T.E.),  perhaps  midchord,  on  the  outside  duct 
surface. 

EFFICIENCY  OF  DUCT  AND  PROPELLER  SYSTEM 

The  theoretical  propulsive  performance  of  the  ducted  propeller  can  be  analyzed  by 
using  the  thrust  and  power  coefficients,  which  were  computed  from  the  propeller  lifting-line 
programs  discussed  previously,  and  estimated  hull  efficiency  factors.  In  terms  of  the  sec¬ 
tion  drag-to-lift  ratio  r  ,  the  hydrodynamic  pitch  angle  /3 f- ,  and  the  nonviscous  propeller  thrust 
and  power  coefficients,  the  viscous  propeller  thrust  and  power  coefficients  are: 

f1  d(CThsh 

CThs=  J  (l-*tan0f) - — - dx  (3) 

*  u 


The  propeller  wake-adapted  propeller  efficiency  r/B  is  the  ratio  of  the  propeller-thrust  power 
(output)  to  the  propeller-absorbed  power  (input)  and  is  given  by 

CTP 

vB  =  7 —  (6) 

LPs 

Finally,  the  total  efficiency  of  the  duct  plus  propeller  is 

/  Td\ 

VT-VB  \l  +  T  J  (7) 

Figure  10  shows  versus  model  rpm  as  obtained  from  Equations  (4),  (5),  and  (6) 
with  the  total  wake  (hull  plus  duct)  as  input  to  the  composite  ducted-propeller  computer  pro¬ 
gram  for  the  given  propeller  diameter  and  thrust-loading  coefficier,'  It  may  be  seen  that  opti¬ 
mum  efficiency  occured  at  800  rpm. 
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The  quasi-propulsive  coefficient  is  defined*  as 


~  p  =  Vt  •1?// 


(8) 


where 

7// 


M 

\ 


(l-<)/(l-»o) 

1-t  is  the  thrust-deduction  coefficient,  and 
l-wQ  is  the  effective  inflow-velocity  ratio. 

From  data  contained  in  Reference  4,  1  -  t  was  estimated  to  be  0.88.  If  it  is  assumed  that 
(1  -  u>0)  -  (1  -  iVy)  where  the  volume-mean  inflow-velocity  ratio  is  defined  by 

r- 

1  V 


i  -  «V  =  ■ 


xdx 


1  -  x. 


(9) 


then  the  hull  efficiency  can  be  estimated  by  an  integration  of  the  curve  of  (1  -  wx)  total 
versus  x  given  in  Figure  6.  A  value  r\^  =  0.88/0.823  =  1.07  was  obtained.  The  theoreti¬ 
cally  predicted  quasi-propulsive  coefficient  for  the  final  ducted  propeller  design  may  be  cal¬ 
culated  from  the  optimum  value  t]b  -  0.675  (Figure  10  and  Equations  (7)  and  (8)).  Thus, 

7jD  =  0.675  (1  +  0.04)  1.07  =  0.75 
where  the  ratio  T d  /  Tp  is  taken  as  0.04. 

PROPULSION  TEST  RESULTS  AND  DISCUSSION 

Figure  11  shows  the  results  of  the  submerged  overload  and  underload  propulsion  test 
of  NSRDC  Model  4620  with  the  duct  and  Propeller  4271.  From  the  experimental  curves  of 
Figure  11,  it  is  possible  to  analyze  and  to  compare  the  powering  performance  of  the  cruising, 
ducted-propeller  system  to  its  theoretically  predicted  performance  and  to  the  performance  of 
an  optimum  wake-adapted,  unducted  stern  propeller. 

Two  generally  recognized  approaches  were  considered  for  interpreting  the  results  of 
the  resistance  and  propulsion  tests.  On  the  one  hand,  the  towed  resistance  could  be  taken 
as  the  bare  hull  and  the  propulsor  could  be  viewed  as  the  complete  duct  and  propeller  unit. 
On  the  other  hand,  the  towed  resistance  could  consist  of  the  hull  plus  duct  without  propeller. 
In  the  latter  case,  the  duct  would  be  viewed  as  an  appendage  to  the  hull.  Both  approaches, 
of  course,  would  lead  to  the  same  prediction  of  shp  and  rpm  for  the  prototype.  However,  in 
the  present  case,  the  author  preferred  to  include  the  duct  in  the  towed  resistance  because  it 
is  easier  to  compare  theory  and  experiment. 


‘The  drag  is  defined  as  the  hull  plus  duct  with  the  various  propulsion  quantities  estimated  as  indicated. 
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Table  3  and  Figure  12  give  the  pertinent  data  for  powering  a  200-ft  prototype  of 
Model  4620  with  a  correlation-allowance  coefficient  ACy  =  0.5  x  10“  3  .  The  data  of  Table  3 
are  given  for  a  constant  ship  speed.  However,  the  various  propulsion  coefficients  are  essen¬ 
tially  constant  over  a  wide  range  of  ship  speed;  V  =  29.3  knots  follows  from  the  design  con¬ 
dition.  It  is  found  from  the  shp  given  in  Table  3  that  the  ducted  propeller  absorbs  about  28 
percent  more  power  than  does  the  unducted  optimum  single-screw  propeller.  The  percentage 
applies  over  the  entire  submerged-speed  range.  The  relatively  low  efficiency  of  the  ducted 
propeller  is  due  in  part  to  its  smaller  diameter  and  higher  rpm,  and  to  the  zero  value  of  Td/T 
at  ship  propulsion. 

TABLE  3 

Comparative  Propulsion  Data  for  a  200-Foot  Prototype 
of  Model  4620  at  V  =  29.3  Knots 

(103ACy=  0.5;  propellers  at  Station  x  -  0.966) 


Quantity 


Estimated  Optimum 
Single-Screw  Propeller* 

(D  =  13.63  ft  Z  =  5) 


Ducted  Propeller 
(Dp=  10.04  ft**  Z  =  7) 


Mp  or  (M)^  0.88  0.88  0.912  0.943 

(1— w«) 

_  0  0.79  0.823  0.823  0.823 

See  Equation  (9) 

♦From  Reference  4. 

♦♦From  Reference  1. 

***77  Cannot  be  compared  for  single  screw  and  ducted  propeller  because  /V-  is  different. 
D  Comparison  should  be  based  on  shp. 


Theoretical 

(Predicted) 

Experimental 

(Actual) 

Experimental  at 
(CThs)  Design 

220  (optimum) 

237 

233 

8900 

8990 

8936 

1.344 

1.246 

1.267 

0.400** 

0.423 

0.400 

0.039 

0.00 

-0.046 

0.70 

0.67 

0.66 

0.675 

0.67 

0.69 

0.75 

0.74 

0.747 

1.07 

1.11 

1.14 

0.88 

0.912 

0.943 

0.823 

0.823 

0.823 

As  mentioned  previously,  the  limited  objective  of  the  present  work  was  to  evaluate  the 
ducted  propeller  configuration  developed  in  Reference  1.  For  the  wake-adapted  case,  the 
best  propeller  is  not  necessarily  that  with  the  largest  practical  diameter.  To  optimize  on  the 
basis  of  minimum  shp,  the  variation  of  the  hull  efficiency  with  propeller  diameter  must  be  con¬ 
sidered.  It  is  very  likely  that  improved  performance  could  be  expected  for  the  present  ducted 
propeller  if  a  larger  diameter  (probably  close  to  that  for  the  optimum  unducted  propeller)  were 
chosen.  An  estimate  indicated  that  the  predicted  4-percent  duct  thrust  would  probably  be  lost 
in  overcoming  the  drag  of  the  four  supporting  struts,  which  was  not  considered  in  the  theory. 

In  attempt  to  further  explain  the  28  percent  more  power  required  by  the  present  ducted 
propeller  than  by  an  unducted  optimum  single-screw  propeller,  propulsive  performance  was 
examined  for  other  apparently  successful  ducted-propeller  configurations.12  It  was  found  that 
the  usual  cruciform-stern-appendage  arrangement  for  submarines  could  increase  hull  efficiency 
by  about  3  percent.  The  configurations  examined  had  utilized  stator  vanes,  and  calculations 
showed  that  their  use  with  the  subject  ducted  propeller  could  improve  total  efficiency  by  more 
than  5  percent.  These  factors,  combined  with  the  stringent  circumstance  of  comparing  with  an 
optimum  unducted  propeller,  go  far  in  explaining  the  relative  propulsive  performance.  It  seems 
likely  from  the  foregoing  that  an  improved  cruising-ducted  propeller  could  be  developed  with 
performance  perhaps  comparable  with  existing  designs. 

It  is  important  to  compare  the  theoretically  predicted  (J  =  1.344)  and  experimental 
(J s  =  1.267)  performances  of  the  ducted  propeller  (Table  3)  at  the  design  thrust-loading  coef¬ 
ficient  C Ths  =  0.400.  In  analyzing  this  result,  it  was  found  that  the  propeller  was  under¬ 
pitched  due  to  an  error  in  input  to  one  of  the  computer  programs  and,  therefore,  that  it  operated 
at  a  rpm  that  was  about  6  percent  high.  Performance  at  ship  propulsion  might  have  been  bet¬ 
ter  had  the  required  thrust  been  known  more  accurately.  In  the  design  process1  model 
boundary-layer  characteristics  were  used,  and  a  ship-model  correlation  allowance  was  not 
considered.  It  is  seen  that  at  design  C  Ths  the  theoretical  and  experimental  values  of  shp 
were  8900  and  8936,  respectively. 

Two  features  of  the  performance  of  the  present  ducted  propeller  at  unusually  high  over¬ 
load  (about  three  times  C t  for  the  ship)  are  worthy  of  comment: 

1.  K  T  increased  almost  linearly  in  the  high  range  of  C , ,  and  the  duct  thrust  was 

1  total 

approximately  40  percent  of  the  propeller  thrust. 

2.  Although  considerable  thrust  was  carried  by  the  duct,  this  portion  of  the  total  thrust 

was  not  delivered  effectively  to  the  hull  as  indicated  by  the  decreasing  total  thrust-deduction 
factor  (l-0totai  •  The  ^act'  ^^ed  duct  augmented  the  hull  drag  relatively  more  than 

did  the  propeller  load  may  be  due  to  the  fact  that  the  center  of  thrust  on  the  duct  shifted  for¬ 
ward  at  large  overloads. 


10 


CONCLUSIONS 


1.  Agreement  was  good  for  the  theoretical  pressure  distribution  on  the  outside  of  the  duct 
as  calculated  by  the  present  method  and  that  of  Kriebel  and  Mendenhall1  for  a  duct  of  zero 
thickness.  Agreement  was  not  as  good  for  the  pressure  distribution  on  the  inside  of  the  duct. 

2.  The  theoretically  predicted  shp  was  within  1  percent  of  the  experimental  shp.  The 
experimental  propeller  rpm  was  6  percent  higher  than  the  design  value  because  the  propeller 
was  underpitched. 

3.  A  200-ft  prototype  equipped  with  the  cruising,  ducted-propeller  unit  required  28  per¬ 
cent  more  shp  than  did  a  conventional  unshrouded  optimum  single-screw  over  the  entire  speed 
range.  The  ducted  propeller  performance  could  be  improved  by  a  parametric  study  to  determine 
the  optimum  propeller  rpm-diameter  relationship  and  other  design  refinements. 

4.  At  high  overload,  the  duct  thrust  was  about  40  percent  of  the  propeller  thrust;  however, 
the  thrust  was  not  effectively  delivered  to  the  hull  as  indicated  by  the  decreasing  total  thrust- 
deduction  factor. 
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Figure  2  -  Comparison  of  Hull  Shapes  for  NSRDC  Model  4620  and  for  the 
Point  Source  and  Line  Sink  of  Reference  l 


Figure  3  -  Cruising  Duct  2.1 
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Figure  6  -  Hull  Wake  and  Total  Wake  of  Hull  Plus  Duct 
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Figure  9  -  Pressure  Distribution  on  Duct 
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Figure  10  -  Ducted  Propeller  Efficiencies  versus  RPM 
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Figure  12  -  Power  and  RPM  Curves  for  a  200-Foot  Prototype  of  Model  4620 
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